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A simple quasiparticle model, motivated by lowest-order perturbative QCD, is 
proposed. It is applied to interpret the lattice QCD equation of state. A reasonable 
reproduction of the lattice data is obtained. In contrast to existing quasiparticle 
models, the present model is formulated in dynamical rather than thermodynamical 
terms, and is easily applicable to a system with finite baryon density. In particular, 
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the model simulates the confinement property. 
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. I. INTRODUCTION 



The most fundamental way to compute properties of strongly interacting matter and, 
in particular, its equation of state (EoS) is provided by lattice QCD calculations 



y, lattice data on the EoS 
4j. Interpretation of these 
is hardly possible in view 



The technique of these calculations rapidly progresses. Recent 
at finite baryon chemical potential became available P, [j" 
data within the straightforward QCD perturbation theory 
of its extremely poor convergence for any temperature of practical interest. To over- 
come this poor convergence, resummation schemes were developed. A scheme based on 
hard-thermal-loop (HTL) effective action [f| have been proposed, with alternative formu- 
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lations in the form of so-called HTL perturbation theory [7[ or based on the ^-derivable 
approximation 8(. This approach justified a picture of weakly interacting quasiparticles, 
as determined by the HTL propagators, and resulted in remarkably good agreement with 
lattice data above 3Tc, with Tq being the critical temperature of the phase transition. It 
is important to emphasize that this quasiparicle picture emerges directly from the QCD 
dynamics, although treated within the thermal framework. Recently, a new resummation 
scheme based on dimensionally reduced screened perturbation theory (DRSPT) was pro- 
posed 0. In certain sense, the efficiency of this DRSPT scheme even surpasses that of 
the HTL perturbation theory. However, it still gives reliable results only above 2.5Tc- 
To extend thi^^vej— below 3T C , various phe_.og.ca. qUaS i- 



particle models (lU llll ll^ llS Il4 U4 lla 

j were proposed. These models are formulated in 
terms of massive quarks and gluons and are constructed in such a way that at high tem- 
peratures they match the perturbative results and then extend them down in temperature. 
It is not clear, if a quasiparicle picture is relevant below 3Tc at all. Therefore, all these 
models are purely phenomenological. Nevertheless, interpretation of lattice data within 
these quasiparticle models turned out to be very successful. With few phenomenological 
parameters it was possible to reasonably reproduce all lattice thermodynamic quantities. 
The feature of all above cited quasiparticle models, which still looks slightly irritating, is 
that they are formulated in terms of thermodynamic quantities (i.e. temperature T and 
baryon chemical potential /x) rather than dynamical ones, like various densities. From 
the theoretical point of view, the quasiparticle picture should be formulated in dynamical 
terms. The thermal equilibrium is only a particular case of this general picture. From the 
practical point of view, if the quasiparticle is introduced as a dynamical object, it would 
be possible to use a quasiparticle model for extending the equilibrium lattice description 
to (at least, slightly) nonequilibrium configurations, relying on reasonable reproduction 
of equilibrium properties by this model. Such kind of extension is really required for 
analyzing heavy-ion collisions, where the thermalization is still a debating problem. 

In the present paper we propose a simple quasiparticle model formulated in dynamical 
terms. In construction of this model, we proceed from properties of the perturbative 
solution to QCD, justifying the quasiparticle picture, rather than from first principles of 
QCD. 
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II. DYNAMICAL QUASIPARTICLES 



Here we will follow the line of Refs. 
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assuming only massive quasiparti- 



cles and avoiding artificial reduction of quark-gluon degrees of freedom as in Let the 
effective Lagrangian for transverse gluons <p a and quarks i[) qc of Nf flavors be as follows 



1 N " 
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0=1 

N c N f 
c=l <j=l 

where iV c = 3 is number of colors, N g = 2(N 2 — 1) is number of transverse gluons, taking 
into account two transverse polarizations, m g and m q are effective masses of gluons and 
quarks, respectively, depending on self-consistent fields 77 and £ = {£1, ...,^ Nf }. B(r},£) is 
a potential of mean-field self-interaction. Writing down Lagrangian {1} we have omitted 
kinetic terms of the rj and £ fields, assuming that they are not essential for the problem. 
Note that these kinetic terms are precisely zero in the spatially homogeneous equilibrium 
and hence are really negligible for slight deviations from it. This Lagrangian is written 
proceeding from general features of the perturbative solution to QCD, which claims that 
quarks and transverse gluons are weakly interacting quasiparticles. Here all interactions 
between gluons and quarks, as well as their self-interactions, are hidden in their effective 
masses depending on mean fields, which in their turn are determined in terms of these 
masses. 

Equations of motion for the mean fields are derived in the standard way: 

,2 N g 1 N f „ 2 7V C 



dri 2 2 dri 2 ^ ™ + 2 ^ dri 2 ^ m Q ' 

£1=1 g=l C=l * 

dB _ Idml^. l^^r^(W 



(2) 



^E«>^ESE™- (3) 



Here 



^M^r « Ha)t ,, (lil)] (5) 
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are scalar densities of gluons and quarks divided by mass, respectively, with f g (x,k), 
f q (x,k) and f q (x,k) being distribution functions of gluons, quarks and antiquarks in 
space-time (x) and 4-momenta (k). In the particular case of thermal equilibrium we 
going to consider here, these are 

fg ^ = expP 2 + m2)i/2/T] -1' (6) 

fgW = exp{[(fc 2 + m2)V2_^]/ T } + l' (?) 
/*(*) ' 



exp{[(A; 2 + m2)i/2 + /ig ]/T} + l- v > 

where T is the temperature, and \i q is the g-quark chemical potential. In general, all \x q 
may be different. If we consider a system with zero overall strangeness and charm, \x q 
relates to the baryon chemical potential as fi u = Lid = yu/3 with all other \x q = 0. 

A solution of Eqs (0) and (J3j), which are usually referred as gap equations, provides us 
with an expressions of fields rj and £ q in terms of above scalar densities. Without loosing 
generality, it is convenient to demand that these solutions for the rj and £ fields are given 
by scalar densities of gluons and quarks divided by mass 

v 2 = E W > (») 

a=l 

^ y{MA_ (io) 

ti rn q 

Indeed, had we started from other collective variables rj and which differ from 77 and 
£ defined by Eqs Q and (JTUJ), and the corresponding potential B(rf,£), equations of 
motion ((2)) and © would provide us with solutions rj(r],£) and £(77, £) with 77 and £ 
associated with densities (jUJ) and IjlOj) . Then we could immediately redefine the potential 
as -8(77, £) = B(rj(r], £), £(77, £)) and thus transform to desired variables 77 and 

The high-temperature limit, T > Tp, where is the temperature of the decon- 
finement phase transition, puts certain constrains on the functional dependence of the 
effective masses m g (r),£) and Tn q (r),£). In this limit, the straightforward calculation of 
above scalar densities in the leading order results in 

r/ 2 (T»T c ) ~ ^T 2 , (11) 
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gfo, T » T c) * N c Qt 2 + ^- 2 4 ) . ( 12 ) 



Perturbative values of m g and m q are also known 



m 2 ({M,T»T c ) = ^ 





/ Nf \ 

(2N c + Nf)T* + ^ti 
\ q=1 ) 



/({/iJ,T»T c ),(13) 



m 2 q ({fx q },T » T ) - m 2 = JJ$L- (V 2 + ^ p 2 ({/x,},T » T c ), (14) 

where m q o is the current mass of the g-quark, Nj G is the effective number of quark flavors 
which can be excited, and g 2 is the QCD running coupling constant squared, generally 
depending on T and all \L q . In the particular case of all fi q = 0, the latter is also known 
\wL \v\ , in the 2- loop approximation it is 

Vi nlT^TM 16tt 2 / 2/3 1 lnln(2 7 rr/A) 2 ^ 



/5 ln(27rT/A) 2 V $ ln(27rT/A) 
with 

A, = ^(UN C - 2Nf), fa = i(34iV 2 - 13iV c iVf + 3Nf/N c ), (16) 

and A being the QCD scale. The energy scale is taken here equal to 27rT, i.e. the first 
nonzero Mastubara frequency. In HTL calculations this scale is sometimes varied from 
ttT to AttT to determine theoretical error bars. 

Expressing T and all fi q in terms of scalar densities rj and £ of Eqs. (jllj) and ()12j) . 
and substituting them into expressions for asymptotic effective masses (fT3)l and (JHJ), we 
arrive at the following expressions for the latter 



"fat) = f^+^E^v^), (17) 

Here we keep Nf instead of Nf, since quark densities £ 2 automatically take care of quark 
contribution, i.e. £ 2 ~ 0, if the g-quark is too heavy. These expressions give at the same 
time the ansatz for the effective masses in terms of scalar densities. 

Now our goal is to find an appropriate expression for g 2 (rj,£), which takes the limit 
(|15|) at fi q = and T ^ Tc, and then to solve gap equations (J2J) and © with respect 
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to B(r),£). Then the quasiparticle model would be completely defined, and we could do 
any calculations both in thermodynamics and nonequilibrium. Unfortunately, we failed to 
solve this problem in general, i.e. we have not found an appropriate potential B which is 
required for calculation of the thermodynamic quantities, cf. Eqs (J27|) and (|2*Hj) . However, 
we have found an elegant solution in the particular case, when number of colors equals 
number of flavors, i.e. Nf = N c . In fact, this case is quite general for comparison to 
lattice data as well as for possible applications in astrophysics and heavy-ion physics. 



III. PARTICULAR CASE OF N f = N c 



Let us consider a particular case when Nf = N c while the quarks may be different, i.e. 
their current masses m q Q as well as chemical potentials [i q may differ. 
Since Nf = N c , Eqs. (JT7J) and (fTHj) can be represented as follows 

Nf 

< = f£c,V(x), (19) 

" q=l 

m 2 q -m% = g^rCyCx), (20) 



' c 



where 



N 



N n 

~l 

1/2 



x 2 = I EC J • (22) 

Here we just guessed that the \ dependence of g 2 is the proper one. This functional 
dependence is required to define the potential B, with which gap equations (J2J) and (J3J) 
give solutions for masses precisely in the form of Eqs (fT9|) and (|20|) . Indeed, gap Eqs. (j2j) 
and (jHI) in terms of new variables ( q read 

^8B_ l dm 2 g 2 1 ^ d{m 2 q - m 2 q0 ) 
dC 2 2dC 2V + 2^ d( 2 * 



n 9 \ * ' d x 2 x 2 U C 'J ^ + W c £ \ 5 ^ 2 + c u x 2 x 
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Nr 



wA 92+x2i ie Ml (23) 



9=1 

2 ~ 



where i = u,d or s. In fact, this is the main trick advanced by the peculiar case of 
Nf = N c . This relation implies that the potential of mean-field self-interaction B(j], £) is 
in fact a function of a single variable Xi an d g a P equations (J2J) and (jSJ) are reduced to the 
single one 

dB 1 2 d( x 2 g 2 ) , nA . 
dx N g dx 

integration of which is straightforward 

B(X) = B c -±- [xY( X ) ~ xWiXc)] + ^ f d XlX i9 2 (Xi) (25) 

iV 3 iV 9 Jxc 

with 5c being an integration constant. In fact, this function B(x) has a meaning of the 
bag constant of the bag model. However, we will refer Be as the "bag parameter" , since 
it is really constant. 

Thus, we succeeded to determine the B potential, which is required for thermodynam- 
ically consistent calculation of thermodynamic quantities (J27|) -()29 |) . We cannot claim 
that this is the only possible solution for this B, because it was obtained as a result of 
certain guess. However, this solution works quite well in reproducing lattice data, as it is 
demonstrated in the next section. 

A reasonable ansatz for the coupling constant itself is as follows 

1571-2 

9 (x) = W?T« W (26) 

where Xc anc ^ Xo are some phenomeno logical parameters, and an auxiliary function f(x), 
meeting the condition /(% — ► oo) — > 1, helps us to choose between 1-loop {fix) = 1) an d 
2-loop asymptotics of the coupling constant, cf. Eq. (|15j). Two reasonable choices of this 
auxiliary function are discussed below, in sect. I1V1 This g 2 (x 2 ) indeed takes the limit 
(fTHJ) at fi q = and T ^> Tc, provided Xo ^ Xc an d properly defined Xc * n terms of A. 
The proper definition in the case, when the temperature is much larger than all current 
quark masses, T ^> m q o, is as follows 



X 2 c 



N Q Nl' 2 ( A ^ 2 
2k 



s 



It is appropriate to mention here that in spite of the declared case Nf = N c , we 
are able to consider less number of flavours, Nf S < 3, cf. Eq. tj!6j). within the same 
formalism. To exclude a g-quark flavor from the treatment at certain temperature T, 
we should take its current mass to be large: m q o ^> T, which implies m q ^> T, cf. Eq. 
(|20|) . In this limit the respective density — ► and simply falls out of the corresponding 
Co, cf. Eq. (121)) . and hence out of the calculation scheme for lighter particles. Taking 
into account that the contribution of this heavy quark into thermodynamic quantities 
(J2ZJ)-(j2IIj) is negligible as compared with that of lighter particles, we see that this heavy 
quark turns out to be completely switched off from the calculation, as if it does not exist. 
If we consider m q0 ^> T ^> Tq, the contribution of this heavy quark disappears even from 
asymptotic formulas fT3j) and ijHJ) . However, the delicate feature of the present solution 
is that we still should keep the flavor summation in Eq. (fTT?|) running though all 3 flavors, 
in order to obtain the proper gluon contribution into the gluon mass, even in the case of 
m qo ^ T ^> Tq. The reason is that for the heavy quark we still have ( 2 = rj 2 , i.e. the 
gluon density which is nonzero. Thus, e.g. for the 2-flavour case, we should keep Nf = 3 
whereas take Nf S = 2 in Eq. (|TB|) required for definition of the coupling constant ([26)1 . 

To summarize, the procedure of solving the model equations is as follows. First, we 
define all the free parameters of the model (xc, Xo, Be), including the auxiliary function 
f{x)- Given the temperature T and the set of chemical potentials fi q , implicit set of 
equations (JH), ©, (fT§j) - ()2"2"j) and (HHJ) should be solved. As a result of this solution, we 
obtain effective quark and gluon masses and the value of x variable, which is required 
for calculation of B(x), cf. Eq. l23j) . Now, when all the quantities are defined, we can 
calculate the energy density e(T, //), pressure P(T, /i), and baryon density nsiT, /i) as 
follows 

n r°° 

e(T,r) = ^§ J k 2 dk (k 2 + mlfl 2 f g (k) 

; AT r°° 

+ E -T / k 2 dk{k 2 + m 2 q Y' 2 [f q (k) + f q (k)} + B( X ), (27) 

q=l 71 ^° 

p(T)/i) = (&/ (k 2 + m 2 g y/ 2f9{k) 

«— L N f°° k 4 dk 

+£^/ WTW jiim + f,m-B( X ), (28) 
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Note that the thermodynamic consistency is automatically fulfilled in this scheme, since 
we proceed from a proper Lagrangian formulation. 

In particular, we would like to mention that the present model simulates the confine- 
ment of quarks and gluons. When temperature and/or chemical potentials decrease, the 
densities, rj 2 and and together with them the variable x drop down. At some value 
of x t ne argument of ln[(x 2 + Xo)/Xc) * n ^q. becomes very close to 1, and hence 
g 2 — > oo. Thus, there are no solutions to the above equations below certain values of 
temperature and chemical potentials. This can be interpreted as a kind of confinement. 



IV. COMPARISON WITH LATTICE DATA 



Our goal is to fit the above described model to the recent (2+1) flavour lattice data 
for nonzero chemical potentials [2J. To be consistent with these lattice data, we accepted 
current quark masses m U Q = m^o = 65 MeV and m s0 = 135 MeV, which were used in these 
lattice calculations. As we have found out, the actual results of our quasiparticle model are 
quite insensitive to variation of m q o from above lattice values to the "physical" ones m u o = 
m^ = 7 MeV and m s0 = 150 MeV. The model also involves several phenomenological 
parameters: the "bag parameter" Be, cf. (12 5 j) . the "QCD scale" xc and an auxiliary 
function f(x), cf. (j27)j). Another parameter xl, as h was expected, should be taken small 
Xo ^ Xc- m f ac t> h shifts the lower limit of integration in the expression for B(x), 
cf. Eq. (125)1 . from the singular point of the coupling constant, cf. Eq. (I2fij) . and hence 
regularizes the calculation of B(x)- Therefore, it is closely related to the "bag parameter" 
Be, which is an integration constant in the same expression. A change of Xo implies the 
corresponding change of Bq. In all further calculations we take x\ = 0.01xc> an d hence 
the below stated values of Be correspond only to this choice. 

An implicit parameter of our model is the critical temperature Tq, i.e. the temperature 
at which the deconfinement phase transition occurs at jj, — 0. We could identify this 
temperature with that of the end point of the solution discussed above. However, this 
end point is numerically determined not quite reliably because of the singular behavior 
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of the solution near it. Another reason is that the end-point temperature should not 
necessary coincide with Tq. The phase transition at fi = in the case of (2+1) flavours 
is of the cross-over type, as it was found in lattice calculations. This implies that a 
strong interplay between quark-gluon and hadronic degrees of freedom occurs near Tc, 
which actually determines the Tc value itself. As we completely disregard the hadronic 
degrees of freedom in the model, we cannot count on proper determinations of Tq value. 
Therefore, we vary Tc from the determined end-point temperature to slightly below in 
order to achieve the best fit of the lattice data. 

As for the auxiliary function /(x), our first choice was 

AW*) = 1, (30) 

which we refer as "1-loop" choice, because with this /i_ioo P the coupling constant takes 
the 1-loop perturbative limit at T — > oo, cf. Eq. (115(1 . In this case we are left with only 
two basic parameters, xc an d B C - These are fitted to reproduce the form of the pressure 
as a function of temperature at zero chemical potential. However, these two parameters 
does not allow us to reproduce the overall normalization of the lattice pressure. With 
this respect, it is suitable to recollect that the overall normalization of the lattice data 
is somewhat uncertain. Indeed, the lattice calculations were done on lattices with N t = 
4 ternpora! extent g To —a the raw lattice data into „ oaes, i.e. to 
extrapolate to the continuum case of N t — > oo, the raw data are multiplied by "the 
dominant T — > oo correction factors between the N t — 4 and continuum 0.518 
and c M = 0.446, [2]. These factors are determined as a ratios of the Stefan-Boltzmann 
pressure at \i = (c p ) and the /i-dependent part of the Stefan-Boltzmann pressure (c M ) 
to the corresponding values on the Nt — 4 lattice [2]. In view of this uncertainty, it is 
legal to apply an additional overall normalization factor to the same quantities calculated 
wi thin _tieie m odel H. m order to k ee P the ournber of ntting par— as few 
as possible, we use a single normalization factor instead of two different ones, c p and c M , 
in the lattice calculations. For the best fit of the lattice data the overall normalization 
factor was chosen equal 0.9 and Tc was shifted slightly below the end-point temperature, 
which by itself was determined quite approximately. Note that the fitted Tc = 195 MeV 
is slightly above its lattice value 175 MeV. The set of parameters is summarized in Table. 
The result of the fit is presented in Fig. ^ In the same figure, also comparison with 
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2-flavour lattice data is presented. For the present "1-loop" variant, 2- flavour data 
are perfectly reproduced with the same set of parameters as for (2+l)-flavour case, only 
the current mass of the strange quark was taken m s0 = 100 GeV in order to suppress its 
contribution. In this case the result for critical temperature is even better: T^ J = 175 

n 

MeV, which well complies with its lattice value |1J. 



Version 


Nf* 


T , MeV 


Xc, MeV 


Bc/Xc** 


/-factor*** 


overall normalization factor 


"1-loop" 


2+1 


195 


141.3 


-97.5 


1 


0.9 


"1-loop" 




175 


141.3 


-97.5 


1 


0.9 


"2-loop" 


2+1 


195 


119.6 


-267.5 


2.6 


1 


"2-loop" 


2**** 


175 


119.6 


-262.0 


2.6 


1 



* Please, do not confuse it with Nf, which should be always Nf = 3. 
** These Be values correspond to the Xo = 0-OlXc choice. 

*** This is the effective value of the auxiliary function f(x) in the temperature 

range under investigation, i.e. from To to 3Tq- 

**** m S Q = 100 GeV in order to suppress the s-quark contribution. 



Table: Best fits of quasiparticle parameters to the lattice data 




Now, when all the parameters are fixed to reproduce the lattice pressure at /x = 0, 
all other calculations can be considered as "predictions" of the model. These results are 
presented in Figs|2HU The model perfectly reproduces /x-dependent part of the pressure, 
AP = P(T,fi) - P(T,fi = 0), Fig. and the "interaction measure", e - 3P, Fig. El 
at various chemical potentials. In particular, it describes practical /x-independence of the 
right slope of e — 3P. At the same time, the lattice baryon density, see Fig. 0J turns 
out to be somewhat overestimated by the model. In fact, this is not surprising, since the 
thermodynamic consistency of continuum lattice limit is somewhat unbalanced because 
of application of different normalization coefficients to the raw lattice data: c p and c M Q| 
mentioned above. Therefore, an exactly thermodynamically consistent model is unable 
to simultaneously reproduce all the continuum lattice data. 

Note that the model fits the lattice quantities only above To- In view of arguments 
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FIG. 1: Pressure normalized by T 4 as a function of T/T c at fi = 0. The solid line corresponds 
to "1-loop" calculation with the overall normalization factor of 0.9. The dashed line represents 
the "2-loop- calculation. The (2 + l)-flavonr lattice data Q are displayed by open circle,, and 
the 2-flavour lattice data |l( - by grey band. 



of Ref. |2Cj this is not surprising. In |2Cj it is argued that below T c these quantities are 
quantitatively well described by the resonance hadronic gas. From this point of view, we 
cannot count on proper description below Tc, since the hadronic degrees of freedom are 
completely disregarded by the model. 

On the other hand, we are able to reproduce the lattice data without varying the 
overall normalization. However, for this we need nontrivial auxiliary function f(x)- One 
of possible choices is 



/woop(x) = 1 + arctan 

with 



-9 (x) In 



A = 0.001- 



16tt 2 



(31) 



which we refer as "2-loop" choice, because with this /2-ioop the coupling constant takes 
the 2-loop perturbative limit at T — > 00, cf. Eq. (I15jl . The additional g 2 ln(0.001) term is 
subleading as compared to g 2 (x) In 9 2 {x) an d thus does not prevent agreement with the 
2-loop approximation for coupling constant. In fact, the function f%i oop is an "exotic" 
representation of a constant, since in the temperature range under consideration, from 
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FIG. 2: AP = P(T, /x) - P(T, fi = 0) normalized by T 4 as a function of T/T c at n = 100, 210, 
330, 410, 530 MeV (from bottom to top) within the "1-loop" (left panel) and "2-loop" (right 
panel) calculations. The (2+l)-flavour lattice data are from Q|. 



Tc to 3Tc, it is /2-i op(x) — 2.6 with good accuracy. Precisely this enhancement of the 
coupling constant is required to fit the actual overall normalization of the lattice data. 
The reason of using function instead of the constant is only that the function provides 
us with the proper 2-loop asymptotic limit. In this respect, any function /, providing us 
with additional factor 2.6 in the temperature range from Tc to 3Tc and respecting the 
proper asymptotic limit of the coupling constant, is as well suitable for this fit. The fitting 
procedure in this case is completely similar to that for the "1-loop" choice. The obtained 
sets of parameters for (2+1)- and 2-flavour cases are also summarized in Table. The result 
of fitting the pressure at fi = is presented in Fig. ^ The 2-flavour case requires here 
only slight tune of the Be value as compared to the (2+l)-flavour case. "Predictions" 
of the "2-loop" version are demonstrated in Figs|2HU The quality of reproduction of the 
lattice data here is approximately the same as in the "1-loop" case. 

In spite of the similar reproduction of lattice data, the two versions of the model reveal 
quite different "internal" quantities, see Figs El and El Their absolute values differ by 
approximately 30%, while the T and \x dependences are very similar in the "1-loop" and 
"2-loop" versions. General trend of these dependences is quite similar to those in the 



thermodynamic quasiparticle model 



14j . Apparently, precisely this trend is essential for 
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FIG. 3: Interaction measure, e — 3P, normalized by T 4 as a function of T/T c at \i = 0, 330, 
530 MeV (which are hardly distinguishable between each other) within the "1-loop" (left panel) 
and "2-loop" (right panel) calculations. The (2+l)-flavour lattice data j^] for different /x = 0, 
330, 530 MeV are displayed by open circles, diamonds and full circles, respectively, similarly to 
other figures. 

reproduction of lattice data within both thermodynamic and present quasiparticle models. 



V. SUMMARY AND OUTLOOK 



We have presented a simple quasiparticle model aimed to interpret the lattice QCD 
data. Similarly to existing quasiparticle approaches , this model is mo- 

tivated by the lowest-order perturbative QCD. However, contrary to those models , it 
is formulated in dynamical rather than thermodynamical terms. Presently we have suc- 
ceeded only for the case Nf < N c , where Nf and N c are numbers of quark flavours and 
colors, respectively. Nevertheless, this is quite a general case for practical applications. 

The model has been applied to fit the lattice (2+l)-flavour QCD EoS at finite baryon 
chemical potentials [2|. This is the most physical and important from the point of view of 
practical applications case. However, we fragmentary considered also 2-flavour case (3]. 
It is demonstrated that a reasonable fit of the quark-gluon sector can be obtained with 
different sets of phenomenological parameters. The "1-loop" version of the model, cf. Eq. 
(HJUJ), certainly looks more natural, since it does not involve an "exotic" auxiliary function 
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FIG. 4: Baryon density normalized by T 3 as a function of T/T c at /x = 100, 210, 330, 410, 
530 MeV (from bottom to top) within the "1-loop" (left panel) and "2-loop" (right panel) 
calculations. The (2+l)-flavour lattice data are from |2j. 
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FIG. 5: Effective masses of « and d quarks (left panel) and s quark (right panel) as functions 
of T/T c at \i = 0, 900, 1500 MeV (from bottom to top). Solid and dashed lines correspond to 
"1-loop" and "2-loop" calculations, respectively. 

/, as it does in the "2-loop" version, cf. Eq. (jUljl . The only problem with the "1-loop" 
version is that it overestimates all lattice quantities by approximately 10% (and slightly 
more for the baryonic density). However, since the overall normalization of the lattice 
data is indeed somewhat uncertain because of the poor extrapolation of these data to the 
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FIG. 6: Effective gluon mass (left panel) and coupling constant squared (right panel) as a 
function of T/T c at [i = 0, 900, 1500 MeV (from bottom to top in left panel, and from top to 
bottom in right panel). Solid and dashed lines correspond to "1-loop" and "2-loop" calculations, 
respectively. 

continuum limit, this misfit is quite acceptable. 

In spite of the difference in absolute values, "internal" quantities of the model, like 
effective quark and gluon masses and coupling constant, reveal very similar behaviour 
as functions of temperature and chemical potential in both "1-loop" and "2-loop" ver- 
sions. Moreover, this behaviour is also similar to that in thermodynamic quasiparticle 



models 



10 



11 



QQ0 



15 



16| ]. Apparently, precisely this general trend is essential 



for reproduction of lattice data within both thermodynamic and present quasiparticle 
models. 

The presented model simulates the confinement of the QCD. In the equilibrium case 
considered here, the solution to the model equations simply does not exist below certain 
combination of the temperature and the chemical potential. In particular, this is the 
reason why we are able to fit the lattice quantities only above Tq- In [20] it is argued that 
below Tc these quantities are quantitatively well described by the resonance hadronic 
gas. From this point of view, we cannot count on proper description below To, since 
the hadronic degrees of freedom are completely disregarded by the model. From both 
theoretical and practical points of view, it is desirable to include hadronic degrees of 
freedom in this model. Then we could count on reproduction of EoS in the whole range 
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of temperatures and chemical potentials. Such a "realistic" EoS would be very useful in 
hydrodynamic simulations of relativistic heavy-ion collisions. 
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